Electron Spin Dephasing due to Hyperfine Interactions with a Nuclear Spin Bath 
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We investigate pure dephasing decoherence (free induction decay and spin echo) of a spin qubit 
interacting with a nuclear spin bath. While for infinite magnetic field B the only decoherence 
mechanism is spectral diffusion due to dipolar flip-flops of nuclear spins, with decreasing B the 
hyperfine-mediated interactions between the nuclear spins become important. We give a theory 
of decoherence due to these interactions which takes advantage of their long-range nature. For a 
thermal uncorrelated bath we show that our theory is applicable down to _B~10 mT, allowing for 
comparison with recent experiments in GaAs quantum dots. 



Introduction. — How a localized electron in a solid 
loses its spin coherence is one of the oldest problems in 
condensed matter physics still attracting attention. It 
has taken on considerable recent significance due to the 
serious efforts in attempting to build a scalable quan- 
tum computer using electron spins in semiconductors as 
qubits [T]. It is therefore extremely important to under- 
stand all aspects of electron spin decoherence in solids 
qualitatively and quantitatively. The hyperfine (hf) in- 
teraction of the electron spin S with a bath of nuclear 
spins Ji (spin-bath) is the main source of the decoher- 
ence of spin qubits at low temperatures, when relax- 
ation effects due to phonons can be neglected. Many 
theoretical studies have been devoted to this problem 
E H 13 EJ H E], and the Spin Echo (SE) decoherence 
time Tse has been recently measured in singlet-triplet [9] 
and single spin qubits [10] based on gated GaAs quantum 
dots. The hf interaction alone can lead to spin decoher- 
ence, either through the direct electron- nuclear (e-n) spin 
flip p], or through the effective (mediated by virtual e-n 
spin flips) intra-bath interaction, which leads to fluctua- 
tions of the Overhauser field felt by the electron through 
the S z Jf part of the hf interaction H0|71|. These 
processes are suppressed at large magnetic fields B, and 
above a certain value of B one has to consider the deco- 
herence due to the intrinsic bath dynamics caused by the 
dipolar interaction between the nuclear spins [3j SJ [5j [6] , 
the so-called spectral diffusion (SD). 

Presently available theories assume the high magnetic 
field limit of £l/A^> 1, where A is the total hf interaction 
energy (sum of the hf couplings with all the nuclei) and 
f2 is the electron spin splitting. This is fulfilled in SE ex- 
periments on electrons bound to the phosphorus donors 
in Si [IT], where the SD theories have shown excellent 
agreement with observations [T]|S][U]. However, in gated 
GaAs dots this limit corresponds to B^> 1 T, while the 
experiments have been done at B 10 — 100 mT and 
therefore the existing theories are of dubious validity in 
this case. The observed 1se< 1 /xs in dots with iV~ 10 6 
nuclei is at least an order of magnitude shorter than the 
theoretical predictions for SD [JJ [5] . This shows that the 



SE decay in the low field regime in GaAs (Sl/A< 1) is not 
due to the spectral diffusion. The existing analytical the- 
ory addressing the SE decoherence due to hf interactions 
only [5] gives negligibly small decay at high B. 

In this Letter we derive a theory of quantum dephasing 
due to the hf interaction, which applies at low magnetic 
fields. Although A/ ft <C 1 is a sufficient condition for 
treating the c-n spin flip as a virtual process |13] , one can 
expect that the necessary condition is weaker. In fact, it 
has been shown 2 that the smallness of the longitudinal 
spin decay is controlled by the parameter S = A/ClyN, 
and the condition of 5 <Cl is much less restrictive, espe- 
cially for dots with many nuclear spins (i.e. large N). We 
argue here that this condition controls the convergence 
of our results at the relevant time scales of interest in 
GaAs spin qubits (as suggested in Ref . [S} , provided that 
the nuclear bath is thermal, uncorrelated, and unpolar- 
ized. Our theory can be used to calculate decoherence 
under any sequence of ideal n pulses driving the qubit 
(e.g. SE corresponds to t/2-ir-t/2 sequence), and we ap- 
ply it to the case of Free Induction Decay (FID) and SE. 
For FID we obtain results which at low fields differ qual- 
itatively from the previously obtained ones [5] [5] • For SE 
we identify the most relevant process contributing to its 
decay: the hetero-nuclear (i.e. involving nuclei of differ- 
ent species) flip-flop of a pair of nuclei. At high B this 
leads to a small oscillation superimposed on the SE sig- 
nal due to the spectral diffusion, while at low B it leads 
to a decay of the SE in ~0.1 — 1/xs in gated GaAs dots. 
This result might explain the recent measurements of the 
SE decay on this time-scale in GaAs [9] 110] . 

Most of the successful theories of spin decoherence 
[H [5] [6] [14] involve an exponential resummation of the 
perturbation series for the decoherence time-evolution 

function W(t) = |p_| (t)\ (where p + _ is the off-diagonal 

element of the electron spin reduced density matrix) . We 
follow an approach along these lines. We take advantage 
of the essentially "infinite range" of the hf-mediated in- 
teractions which couple effectively all the N ~ 10 5 — 10 6 
nuclear spins in the bath (i.e. in the dot) . We consider the 
ring diagrams, identified as the most important contribu- 
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tions for long-range interactions [B] in the perturbation 
expansion of W{i). Corrections to these are suppressed 
as X/N. This approach is closely related to the l/z (z 
being the coordination number) expansion of the parti- 
tion function of the high-density Ising model [TH] , as well 
as recent calculations of the influence of the hf interac- 
tions on the electric dipole spin resonance in quantum 
dots [16]. Taking into account the ring diagrams only, 
we can sum all the terms in the cumulant expansion of 
W(t). 

Effective Hamiltonian. — The original Hamiltonian 
consists of Zeeman, hf, and dipolar interactions. The 
latter will be neglected below, since decoherence due 
to them can be calculated using well-tested methods 
[3]. The Zeeman energies are Hz = QS Z + X^aH^f 
and hf interactions are given by H^f = A;S ■ J; = 
AiS z Jf + V s {, where i labels the nuclear sites and a[i] 
labels the nuclear species at i-th site (assigned randomly 
to sites in case of multiple isotopes, e.g. 69 Ga and 71 Ga). 
V s f ~ J T is the spin-flip part of the hf interaction. 
The hf coupling Ai is proportional to its species' total 
hf energy A a [i\ and the square of the envelope function 
/ i = l*(r l )| 2 : A i = A ot [ i ]f i with normalization Y.ifi = n c 
(n c being the number of nuclei in the unit cell, equal to 2 
in GaAs). The number N of nuclei interacting apprecia- 
bly with the electron is defined by N = J2i fi/Tli fi- For 
any a, J2iea — n aA 2 i /N where n a is the number of a 
nuclei per unit cell, and the maximal A iea mA a /N. In 
calculations below we use a 2D Gaussian wave-function 
(with the parameters for GaAs given in Fig. [lj, but the 
results on time-scales relevant for decay at low fields 
(t <C N/A) are determined only by N, not by detailed 
shape of ^(r). 

We derive the effective pure dephasing Hamiltonian 
H by removing the V^f interaction with a canonical 
transformation H] ■ In the second order with respect 
to V s { we obtain the two-spin (2s) interaction consid- 



ered in Refs. [3 [8]: H, 
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ZS'YntjBijJtJl with 



Bi j = AiAj/AQ. We will refer to this term as the S z - 
conditioncd term. These non-local hf-mediated flip-flop 
interactions are the most important at moderate fields 
and yield the effective Hamiltonian 
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The electron Zeeman term was discarded as being a con- 
stant of motion; the nuclear Zeeman terms, however, 
are important because B it j can couple nuclei of differ- 
ing species (this is crucial for SE). At higher orders we 
will have additional ^-independent and multi-spin in- 
teractions which we address near the end of this Letter. 

Decoherence function. — We write H = Hq+2S z Vi+V2, 
with the flip-flop interactions divided into S^-conditioned 
V\ and 5 z -independent V2. We define the function f(t; t') 
which encodes the sequence of n pulses applied to the 



qubit 17], f( t ;t>)=Y:i= 1 (-i) n+k ®(t' -tk)e(t k+1 -t'), 

where tj- with k=l, n are the times at which the pulses 
are applied, and to = 0, t„+i = t, with t being the total 
evolution time. The decoherence function W(t) can be 
written introducing the Keldysh contour [BJ HU [TH] 



W(t) = 





'-'J 




Jc 



[cf(t; t)V 1 {t c )+V 2 {t c )]o17 



(2) 

where Tq denotes the contour-ordering of operators, (...) 
is the ensemble average, t c = (t, c) is the time variable on 
the contour, with c= + (— ) on the upper (lower) branch 
of the contour, and V\^{t c ) are the interactions in which 
the spin operators are given by 

J j ± ( Tc ) = jfe ±iw i T±ic f° = Jfdf{r c ) . (3) 

One remark is in order for the case of FID, in which the 
above formulas refer to the single-spin decoherence, de- 
fined in [5] where it was shown that in unpolarized ther- 
mal bath one can factor out the inhomogeneous broaden- 
ing from W(f) (this is equivalent to FID calculation using 
a "narrowed" state with no net nuclear polarization [8 ). 

Ring diagrams. — We concentrate now on interactions 
from Eq. (JT|). Perturbation expansion of W(t) is W = 
1 + J2 k =i W^ 2k \ with WW - (V 2k ). We use here the 
thermal nuclear bath at temperature very high compared 
to the nuclear Zeeman energies, so that (...) ^Tr{...}, and 
inside the averages all the J 8 + operators have to be paired 
up with Jf operators. Thus, W^ 2k ^ has the following 
structure: 



(4) 



with at most 2k indices different. The full diagrammatic 
expansion is very cumbersome, since the commutation 
rule [J 4 + , J J] =25ijJ z does not allow the standard form 
of Wick's theorem to work [SI Q3]. However, for "infi- 
nite" range interaction (all N spins coupled comparably 
with each other) we have ~ N 2k terms in Eq. Q . From 
the sum over 2k indices we now take only the terms in 
which none of the indices is repeated. This amounts to 
introducing a 1/N error. The simplification is tremen- 
dous: now we only have one pair of J^f and for given 
i, and they effectively commute, as the averages contain- 
ing Jf as the sole operator of the z-th nuclei are zero 
in uncorrelated and unpolarized bath. We can then get 
rid of time-contour ordering, and classify all the terms in 
the perturbation series in terms of ring diagrams, which 
involve clusters of nuclei connected by two-spin interac- 
tions in a cyclical manner. Introducing the matrix Tij 
we can write the ring diagram containing 2k nuclei as 



R.2k(t) — 2_j -^1*2 



(t)...T i2kil (t)~TrT(t) 



2 k 



(•5) 



where T l0 = (l-Si^^aiajBij J c cf(t;r)df (T c )d 3 (r c )dr c , 
with Oi=I Ji( Ji + 1). With the same 1/N error, we may 
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approximate all of the terms in the expansion of W(t) 
as products of R 2 k defined in Eq. ([5]), and In W(t) is the 
sum of the linked terms [19] , leading to 



N = 10 , g = -0.44 



W(t)=exp[j2 

fe=i 



(-l) fc 
2k 



R2k(t) 



(6) 



Using Eqs. and ([6| we arrive at formula involv- 
ing the eigenvalues A; of the T-matrix: W(t) = JL(1 + 
A 2 ) -1 / 2 . Since depends only on and Ai, we may 
approximate R 2 k by considering continuous distribution 
p(A) of Ai and then coarse-graining it. With the relevant 
ranges of Ai divided into Ma slices and with Nj nuclear 
species, we have to deal with matrices of M x M size with 
M = MaNj. The calculation of W(t) quickly converges 
as we increase M with the coarse-graining error negligi- 
ble for t<^NM A /A. In fact, we will show that M = Nj 
(equal to 3 in case of GaAs) gives a good approximation 
at experimentally relevant short times. 

Free Induction Decay. — The main source of the FID 
decoherence is the interaction from Eq. (JTJ with only 
homo-nuclear flip-flops taken into account (as long as 
Lo a f3 = uj a — ujf3^> A/N, which is fulfilled in GaAs with 
7V«10 6 and BwO.l T). Then we can consider only ring 
diagrams involving only the nuclei of the same species, 
so that Ta 



y ■- sva.{A i jt/2)/A i j, with A^j = Ai — Aj, and 
R2k = Y, a R $k, so that ^FiD(t)=n Q ^ ID (t). For short 
times, when t -C N/A (i.e. t -C 20 fjs for GaAs with 
N = 10 6 ), we obtain a very simple result: R 2k = (R 2 ) fe 
and R 2 = £,at 2 , leading to the resummation of the expo- 
nentiated series in Eq. © giving us 



V^+^t 2 



5« 



n 2 A 4 



(7) 



so that the characteristic decay time is Tfid ~ Nil /A 2 as 
in Ref. \S\ In Fig. [I] we show that this formula describes 
very well the FID decoherence in GaAs for B < 1 T. Let 
us also note that for InGaAs with N= 10 5 , <j e ff = 0.5, and 
at 5 = 6 T from Eq. @ we obtain an estimate TpiD~0.6 
/is which is close to the recently observed value of 3 /is 
in an optically selected "narrowed" state [2U]. Eq. Q 
should be compared with the "pair approximation" result 
of Ref. [5j which in our notation reads W a = exp(— hR% ). 
This agrees with Eq. Q for i? 2 <l, when both formu- 
las give W~ 1 — 2^2; but at times at which R 2 becomes 
comparable to 1 only Eq. ^ correctly takes into account 
the correlations of groups of spins larger than pairs. In- 
terestingly, for long times (t^>N/A), the ring diagrams 
with more than two spins are suppressed by oscillations 
in , and the only remaining term is then R 2 ~t at this 
time-scale. Thus we recover the long-time Markovian be- 
havior discussed previously [SJ [5] . 

Spin Echo. — If we use only the S^-conditioned interac- 
tions and allow only for homo- nuclear flip-flops, there is 
no SE decoherence [5], i.e. WsE(t) = l- The SE decay can 
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FIG. 1: W(t) due to hf-mediated interactions for FID calcu- 
lated for two GaAs dots of different sizes and various B. The 
solid lines are the results of the full T-matrix calculations 
(with convergence achieved at matrix size M«100), dashed 
lines are obtained with 3x3 T-matrix taking into account 
hetero-nuclear interactions, the dotted lines use Eq. 0. For 
GaAs we have used n a =0.604, 0.396, 1, A a = 5.47, 6.99, 6.53 • 
10 10 s _1 , and uj a /B = -6.42,-8.16,-4.58 • 10 7 s^T" 1 for 
nuclei of 69 Ga, 71 Ga, and 75 As, respectively. The electron 
spin splitting is Q = 0.88 g eS B[T] ■ 10 11 s" 1 with cf cff = -0.44 
in a large quantum dot. At high B, the actual decay will be 
bounded by ~10 [is due to spectral diffusion [5]. 



come then either from .^-independent interactions, or 
from considering the hetero-nuclear flip-flops in Eq. (JTJ). 
The latter lead, for t<S^N/A and u) a p ^>A/N, to 



T, 



SE 



A k Ai 



8ie iu af3 t/2 

w a 



sm 



(8) 

The multiplicative dependence of the above matrix on 
Ak and Ai allows us to define a matrix T a p of Nj x 
Nj dimension, in which A k is replaced by y/n a /NA a [k]- 
Then we have R2k(t) = J2i m terms of Nj eigenvalues 
of T a p. In GaAs we have Nj = 3, and we get A, =0, ±A, 
so that W SE (t) = [1 + \R 2 {t)]-\ where R 2 = 2(\f 12 \ 2 + 



\Tiz\ 2 



|T 23 ' 2 



) and 



- 2 4AaA 2 p , 4 LO a[j t 

' /3 ' = N 2 fl 2 LU 2 "" n /3 a " Q /3 Sln A 



(9) 



a/3 



The SE decoherence is shown in Fig. [2] At high B when 
|T ai 3| 2 « 1 we obtain a small oscillation which in ex- 
periments will be superimposed on the SD decay. At 
low B, which for N— 10 5 (10 6 ) in GaAs corresponds to 
n < n c ~ ^/r~A/VN ps 250(75) mT (where r ~ 10 3 is 
the ratio of electron g-factor to the typical difference of 
nuclear (/-factors), the coherence decays to zero in char- 
acteristic time TsE~r/fl c ps 0.25(0.75) /is. This is close 
to T SE « 0.3 - 0.4 fjs observed in 10 at Bk 50 - 70 mT. 
The B dependence of decay is visible in a narrow range 
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FIG. 2: Spin echo decoherence W SE (t) in GaAs due to hetero- 
nuclear processes. The dots are obtained in the A/N <C 
uj af3 , 1/t limit, when W(t) = [1 + ^R 2 (t)]~ 1 , while the solid 
lines are the results of the calculation with the full T-matrix. 
The differences between the two approaches on a microsecond 
time-scale are visible for the smaller dot (upper panel), but 
are negligible for the larger one (lower panel) . The decay time 
due to spectral diffusion (calculated as in Ref. 4) is > 10 us 
in both cases. The GaAs parameters used are the same as in 
Fig.0 



of fields when the decay is substantial but not complete. 
The S z -independent two-spin term. — The contribu- 



tion of H. 



(3) 
2s ' 



'Eij CijJTJJ with a, ocA^Aj/n 2 to FID 



and SE is calculated in the analogous way, and it turns 
out to be negligible for B> 10 mT and iV> 10 5 . 

Influence of multi-spin interactions. — The most im- 
portant multi-spin terms are the ones involving n nu- 
clei in the n-th order in expansion of H, e.g. H. 



(3) 

3s 



S z Y,nk C ijkJ?JiJk with djk oc (A/N) 3 /ft 2 . Apply- 
ing the operator norm \\A\\ = max/$i$\ = i | (<&| A |<E>) | to 
these terms, one gets ||ff (n) || ~ A{A/^l) n - 1 , so that it 
may seem that the expansion is convergent only when 
A/Cl < 1. However, the state |$) maximizing the norm 
is a highly entangled state of nuclear spins. Here we 
are interested in the evolution of the nuclear spin sys- 
tem averaged over an uncorrelated thermal ensemble, and 
therefore such massively entangled nuclear states do not 
appear in the calculations. 

We calculate the lowest order term in the expansion of 

(2) 

W3 S (t) w 1 + W$ s + ... in powers of three-spin (3s) inter- 
action. For both FID and SE we get that for t < N/A 
(when W(t) decays practically to zero due to the two-spin 
interactions) we have /Rq ~ S 2 . Influence of each 
higher-order multispin term in H comes with a higher 
power of the small parameter 6, and our neglect of these 
terms is a well-controlled approximation at the time-scale 
of interest. We therefore believe that the approach used 
in this work remains valid at relevant times down to 
B w 10 mT for GaAs dots with iVwlO 5 - 10 6 nuclei. 



Conclusions. — We have presented a theory of quan- 
tum decoherence due to hyperfine interaction of a central 
spin with a thermal nuclear bath. For both FID and SE 
we obtain analytical formulas for non-exponential time 
dependence of decoherence. For SE we obtain a good 
agreement with the measurements in gated GaAs dots, 
and predict the appearance of characteristic oscillations 
in the signal at magnetic fields slightly higher than the 
ones used in experiments. 
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